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RESUMEN 

En este articulo se procede a mostrar y comparar formas sencillas y complejas de modelos 
de momento curvatura y matrices de rigidez de elementos con secciones de diferentes rigideces en 
un rango lineal o no lineal, enfocado en la programacion para analisis no lineales. 
Al tener los tipos de modelos con los que se puede trabajar se concluye con la comparacion de 
estos y se expone las ventajas y desventajas de cada modelo y tipo de matriz de rigidez para 
analisis no lineales. 



ABSTRACT 

In this paper we proceed to show and compare simple and complex forms of models of 
moment curvature and stiffness matrices when elements have sections with different stiffness in a 
linear or non-linear range, focusing on programming for nonlinear analysis. 
When we have all types of models that we can work, concludes with a comparison of all their 
advantages and disadvantages, for each model and type of stiffness matrix for nonlinear analysis is 
described. 
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1. 



INTRODUCTION 



El diagrama momento curvatura nos permiten conocer la ductilidad de una seccion, es 
decir su potencial de incursionar en un rango no lineal antes de llegar a su fallo, esto se logra a 
traves de conocer el comportamiento del modulo de elasticidad (E) y la inercia (I) (<t> (curvatura) = 
M(Momento)/EI), suponiendo que el fallo sea causado por flexion y no por otro tipo de esfuerzo 
como por cortante. Los codigos sismoresistentes guian a un tipo de fallo por flexion. 



2. COMPORTAMIENTO DEL ACERO ESTRUCTURAL ESFUERZO VS 
DEFORMACION 



En el diseno y construccion de estructuras como edificios y puentes se trabaja con 
materiales de comportamiento quimico y mecanico predecibles limitado por normas locales, tanto 
en composicion como en resultado. El acero estructural tendra limites de: fluencia, maximo 
esfuerzo y rotura dependiendo de la composicion quimica, se puede trabajar con distintos tipos de 
modelos Esfuerzo-Deformacion Unitaria. 

Un modelo muy completo y utilizado por programas de CSI Berkeley para graficar 
diagramas momento curvatura, es el de Holzer de 1975. 



Los siguientes parametros definen el modelo de Esfuerzo Deformacion del acero: 

e=Deformacion unitaria del acero, f=Esfuerzo del acero, E= Modulo de elasticidad, fy= Esfuerzo de 
fluencia del acero, fu= Maximo esfuerzo del acero, esh= Deformacion unitaria en comienzo de zona 
de endurecimiento, eu= Deformacion unitaria correspondiente al maximo esfuerzo del acero, er= 
Deformacion unitaria en la rotura del acero, ey= fy/E. 

Para la zona elastica (e^ey). 




Deformacion Unitaria, £ 



Figura 1 Zonas de curvas Esfuerzo-Deformacion del acero 



f=E*£ 



(1) 



Para la zona perfectamente plastica (£y<£^£sh). 



f=fy 



(2) 
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Para la zona de endurecimiento y decaimiento (£sh<£<£r). 



/ = /y (i + r (g-i) e a-)) 



(3) 



Donde, 



r — 



£ ~ £ sh 



(4) 



£ u £ sh 



Todos los parametros se los puede conocer a traves normas locales que delimitan a estos, 
por ejemplo, muchos edificios son disenados con acero ASTM A36, entonces de la norma ASTM 
A36 podemos conseguir: 

fy= 2530 kg/cm2 
fu= 4080 kg/cm2 
£r= 20% 

Del codigo AISC 360 se puede obtener el modulo de elasticidad: 

E= 2043000 kg/cm2 

Entonces, £y= fy/E= 0.0012=0.12% 

£sh y £u se los puede obtener de ensayos en acero A36, datos que se los pueden comparar con 
bibliografias como del ASM International - The Materials Information Society el libro Atlas of Stress 
- Strain Curves. De aqui se obtuvo: 

£sh= 2% 

£u= 14% 



3. MODELO MOMENTO CURVATURA COMPLETO INCLUYENDO FUERZA 
AXIAL 



Para una seccion sometida a flexocompresion se deben superponer las deformaciones 
debido a carga axial y a momentos. Teniendo estas deformaciones resultantes se procede a 
determinar en que rango, lineal o no lineal, esta cada fibra y asi se obtendran esfuerzos 
correspondientes segun la relacion esfuerzo deformacion del material, este proceso se puede 
observar abajo para varios casos: 
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Para esfuerzo f<=fy: 



X 



Y 



Y 



.XX 




M — Fc* brazo s + Ft * brazo_i 



<D = 



M = Sx * / 

e e 



yi/2 ys/2 



Para esfuerzo f>fsh: 




M — Fc* brazo s + Ft * brazo_i 

£ £ 
<I> = = 

yi/2 ys/2 

Tomando el edificio en dos dimensiones de la figura 2 se obtuvo la Curva de Capacidad 
utilizando para cada seccion una curva de modelo completo a flexocompresion para columnas y 
solo a flexion para vigas. A continuacion se presenta los resultados comparandolos con curvas de 
capacidad sin tomar en cuenta la fuerza axial en el analisis pushover. 
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Figura 2 Estrcutura a ser analizada 

Curva V-Deformaci6n tope 




0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Figura 3 Resultados de curva de capacidad desde diagramas Momento-Curvatura 



NORMATIVA PARA DIAGRAMAS DE MOMENTO CURVATURA 



Para el calculo de analisis no lineales el codigo ASCE 41- Seismic Evaluation and Retrofit 
of Existing Buildings, no da parametros para definir diagramas Momento - Curvatura sino 
diagramas Momento - Rotacion basados en caracteristicas, comportamiento y utilizacion de las 
secciones. 



104 



David Mora, Roberto Aguiar 





t 

b 


5 




. a *| 


/ ° E 


8 or A 



Figura 4 Curva generalizada Fuerza-Deformacion de componentes o elementos 

Para poder trabajar con una relacion EI = M/O como se ha visto se puede proceder con el 
diagrama Momento Rotacion trabajandolo por el metodo Bilineal de Interaccion Carga axial - 
Momento, relacion de momento constante, segmento finito - resorte finito, o momento curvatura - 
segmento finito, etc, tambien se puede utilizar este diagrama de forma directa con las rotaciones 
usando otros modelos. En este articulo se expondra para su utilizacion el metodo Bilineal de 
Interaccion Carga axial - Momento, Modelo de Giberson, y de Matriz de rigidez Condensada en 
coordenadas reducidas, tambien se comparara una aproximacion para relacionar la rotacion con la 
curvatura en el metodo de la longitud de la rotula plastica. 



5. COMPARACION 0y CON EXPRESION ASCE41 

El codigo ASCE 41 da la opcion de usar la expresion general: 

6y = (5) 

J 6*E*I 

Donde, 

E = (modulo de Young) 
I = Momento de inercia. 
L = Largo de elemento 

Z = Modulo plastico, Z = J y dA — Ac * brazojs + At * brazoj 
fy = Esfuerzo de fluencia del acero. 
P = Fuerza Axial aplicada. 
Py=A*fy. 

A = Area de la seccion. 

Esta expresion se obtiene considerando que los dos extremos de un elemento estan 
fluyendo: M=Mp (Momento Plastico en nudo inicial); M'=M'p (Momento Plastico en nudo final). 




Figura 5 Viga deformada por flexion 
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y Ay 



0Ay ~ 6EI 



MpL M'pL 
3EI ~ 6EI 

SiMp=Mp': 

MpL 



'By 



Donde, Mp es el momento plastico. 

Si se toma esta consideracion se obtienen resultados conservadores que pueden ser no 
precisos como se muestra en la siguiente comparacion para el edificio de la figura 2. 

V-Desplazamiento en tope 



200 



50 





































































































Con 


Modelo teta de ASCE 41 

ulando teta y en cada ciclo de carga 










Calc 





















0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 



Figura 6 Comparacion de curvas de capacidad utilizando diferentes formas de caculo 

para 8y 

LONGITUD ROTULA PLASTICA 




Figura 7 Diagrama de momentos de viga de longitud L. 

Para calcular la longitud de la rotula plastica: en la figura 7 se puede ver que los 
momentos maximos seran Mi y Mj, el limite donde fluye la seccion es My, ademas se supone que 
el diagrama de momento es lineal (lo cual puede no ser real si la carga gravitacional es muy 
grande). 

Por triangulos rectangulos podemos determinar Lp: 



Lpa= |M«L L;Lpfc = |Mzi^|, L 



Mi-Mj 



Mi-Mj 



(6) 
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7. MODELOS PARA ELEMENTOS EQUIVALENTES CON SECCIONES DE 
RIGIDEZ VARIABLE 

7.1 Metodo bilineal de interaccion carga axial- momento 

En este metodo la curva de momento curvatura es integrada por dos componentes 
imaginarios, como se muestra en la figura 9. 




Figura 8 (a) Modelo Momento Curvatura Bilineal; (b) Componente Lineal y Elastoplastico 

En esta figura, las pendientes de los componentes lineal y elastoplastico son a1=p*EI, 
a2=q*EI, y p+q=1, donde p es la fraccion de la rigidez flexural de la componente lineal y q es la 
fraccion de la rigidez flexural de la componente elastoplastica. La pendiente de post - fluencia de la 
componente elastoplastica es igual a cero. 




Figura 9 Miembro No Lineal 



Un miembro no lineal tiene dos componentes, el lineal y el elastoplastico. 8i y 8j son las 
rotaciones totales en los extremos del miembro; ai y aj son la rotacion plastica en cada extremo del 
componente elastoplastico. La matriz de rigidez en cualquier ciclo incremental puede ser formulada 
de acuerdo al estado de fluencia. El estado de fluencia puede ser una de las siguientes cuatro 
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condiciones: (a) ambos extremos lineales, (b) extremo i no lineal y extremo j lineal, (c) extremo i 
lineal y extremo j no lineal, y (d) ambos extremos no lineales. 

La matriz de rigidez para la condicion (a) es: 

p-g ap 

La matriz de rigidez para la condicion (b) es: 



fAMi) 



pa pb 1 f A0i 



} (8) 



lAMjj [pb pa + qej {a0j 

La matriz de rigidez para la condicion (c) es: 

[AMi] _ pa + qe pb] f A0i) 

lAMjj - L pb paj lABjJ (y) 



La matriz de rigidez para la condicion (d) es: 

[AMi] 
lAMjj 

Donde: 

a = 4EI/L 
b = 2EI/L 
e = 3EI/L 



Como se muestra en la figura 9, la capacidad del momento nominal, Mn, es usado para 
conocer la condicion de fluencia en el extremo del miembro. Durante el analisis pushover 
monotonico, la capacidad del momento nominal es influenciada por la magnitud de la carga axial. 
Esta simple demostracion nos limita a que todas las secciones a lo largo del elemento deben ser 
iguales- 



7.2 Modelo de Giberson 

En el modelo de Giberson de 1969 se tiene un elemento con tres diferentes rigideces a 
flexion Ela, Elo y Elb y una rigidez axial igual a EA como se muestra a continuacion: 



<>(EI) a (El) 0 (EI) bM 




L 



Figura 10 Modelo de Plasticidad concentrada de Giberson 
La matriz de flexibilidad para este modelo es: 
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F = 



6EI0 



(1 + Sa) -1 
-1 (1 + Sb) 



0 
EIo 

Sa = — ; Sb 
Ela 



0 

EIo 
" Elb 



0 
0 

6EI0 



EA 



Se obtiene la matriz de rigidez como la inversa de la de Flexibilidad: 

Kl = F" 1 



Kl = [(6 * Ela * (Elb + EIo))/(lon * (Ela + Elb + EIo)), 
(6 * Ela * EIb)/(lon * (Ela + Elb + EIo)), 
0, 



(6 * Ela * EIb)/(lon * (Ela + Elb + EIo)), 0 
(6 * Elb * (Ela + EIo))/(lon * (Ela + Elb + EIo)), 0 

0, EA/lon] 

(11) 



7.3 Matriz de rigidez condensada en coordenadas reducidas 

Si no suponemos que la plasticidad estara concentrada en puntos especfficos del 
elemento si no en longitudes Lp (Longitud de la rotula plastica) se puede determinar la matriz de 
rigidez basica de un elemento K1 en un sistema de coordenadas reducido, de la siguiente manera: 




Figura 11 Elemento en coordenadas de sistema reducido 

i) Trabajamos solo en coordenadas reducidas x - y, con tres elementos Lpa, Lo y Lpb de 
diferentes rigideces a flexion Ela, EIo y Eib, y diferentes rigideces axiales EAa, EAo y EAb. 
En la matriz VC se tiene los grados de libertad por elemento. 

L — [Lpa Lo Lpb]; 
EI = [Ela EIo Elb]; 
EA = [EAa EAo EAb]; 

VC = [0 0 1 4 5 6 

45 6789 

789 3 0 2]; 
ngl = max{max(y r C)); 

sen = 0; cosen = 1; 
mbr — 3; 
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ii) Se obtiene la matriz de rigidez en sistema de coordenadas reducidas K1 de cada 
elemento, se le multiplica por matrices de transformacion T12 para llevar al elemento a 
coordenadas locales del elemento K2=T12 T *K1*T12 y a esta por las matrices T23 para 
llevar al elemento a coordenadas generalizadas K3=T23 T *K2*T23. Luego se ensambla la 
matriz de rigidez de la estructura con la matriz de colocacion VC. A continuacion se 
muestra un ejemplo del proceso como un bucle programado en Matlab. 

for i=l :mbr 
lon=L ( i ) ; 



Kl= [4*EI (i) /Ion 
2*EI (i) /Ion 
0 

T12=[0 1/lon 
0 1/lon 
-1 0 



2*EI(i)/lon 0 
4*EI(i)/lon 0 

0 EA(i)/lon]; 



-1/lon 
-1/lon 

0 



0 
1 

0] ; 



T12t= 



[ 0 
1/ (Ion) 

1 
0 

-1/ (Ion) 

0 



0 

1/ (Ion) 

0 
0 

-1/ (Ion) 

1 



-1 

0 
0 

1 

0 

0] 



K2=T12t*Kl *T12 ; % Matriz de rigidez en coordenadas locales 



T23= 



sen 
cosen 
0 
0 
0 
0 



0 
0 
0 

cosen 
-sen 

0 



0 
0 
0 
sen 
cosen 
0 



0 
0 
0 
0 
0 

i] ; 



[cosen 
-sen 

0 
0 
0 
0 

K3=T23 ' *K2*T23; 

% Ensamble de la matriz de rigidez con la matriz de colocacion VC 
for j =1 : 6 

J=VC(i, j) ; 

if J==0 

continue 

end 

for m=l : 6 

mm=VC ( i , m) ; 
if mm==0 

continue 

end 

if J==l && mm==l 

SS ( J, mm) =K3 ( j , m) ; %Matriz de rigidez de los 3 elementos 
SS (ngl, ngl) =0; 

else 

SS ( J, mm) =SS (J,mm) +K3 ( j ,m) ; 

end 

end 

end 



end 



iii) 



Si trabajamos con la matriz total del elemento SS: 
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Q = SS *q, donde, Q es el vector de cargas en cada grado de libertad. 

q = SS' 1 * Q, donde, q es el vector de deformaciones en cada grado de libertad. 



fql 
q2 
q3 
q4 
qS 
q6 
ql 
q8 

Vq9) 



Sinvll 
Sinv91 



Sinvl9 
Sinv99 



Q2 

Ql 
o 
o 

0 
0 
0 

^ 0 J 



Creamos la matriz Sinv=SS _1 y escogemos las tres primeras filas y columnas puesto que 
el resto de cargas del vector Q son cero, ademas estas filas y columnas escogidas son las que 
corresponden a los grados de libertad que necesitamos. Luego volvemos a sacar su inversa 
teniendo: 

Sinvll SinvYl Sinvl3Y~ r 
Kl = Sinv21 Sinv22 Sinv23 
Sinv31 Sinv32 Sinv33. 



ii) 



Se puede obtener el mismo resultado subdividiendo a la matriz SS en una submatriz 
Kn, del tamano de los grados de libertad que necesitamos y procediendo de la 
siguiente manera: 



„„ _ [Kn Bn\ 
~ VBnt Ke\ 

I 0 J Writ Kei U^i 



0 J Writ Kei iqei 
Despejando y simplificando Qn: 

Qn = [Kn — Bn* Ke -1 * Bnt] * qn 
Por lo tanto la matriz en coordenadas simplificadas es: 

Kl = Kn-Bn*Ke~ 1 *Bn 

iv) El resultado que se puede programar sera: 
EC1 

_2EIaEIbEIo(EC5 + ECU + EC10 + EC8 + 6 Ela Elb Lo Lpa Lpb + EC9 + EC3 + EC6 + EC7 + EC4) 
~ EC2 

4 Ela Elb EIo (ECS + ECU + 3EIaEIoLo 2 Lpb + 3 Ela Elb Lo Lpa 2 
_ +6 Ela EIo Lo Lpa Lpb + EC9 + EC3 + 3 Ela EIo Lpa 2 Lpb + EC7 + EC4) 
~ EC2 

4 Ela Elb EIo (ECS + EC1Q + 3 Elb EIo Lo 2 Lpa + 3 Ela Elb Lo Lpb 2 



EC13 = 



+ 6 Elb EIo Lo Lpa Lpb + EC9 + EC3 +3 Elb EIo Lpb 2 Lpa + EC8 + EC 3 + EC6 + EC4) 

EC2 
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EC1A ■■ 



EAa EAb EAo 



EAaEAbLo + EAa EAo Lpb + EAb EAo Lpa 



Kl = 



ECYl EC1 0 
EC1 EC13 0 
0 0 EC14 



(12) 



Donde: 

EC2 = Ela 2 Elb 2 Lo 4 + 4 Ela 2 Elb EIo Lo 3 Lpb + 6 Ela 2 Elb EIo Lo 2 Lpb 2 + 4 Ela 2 Elb EIo Lo Lpb 3 

+ Ela 2 EIo 2 Lpb 4 + 4 Ela Elb 2 EIo Lo 3 Lpa + 6 Ela Elb 2 EIo Lo 3 Lpa 2 

+ 4 Ela Elb 2 EIo Lo Lpa 2 + 12 Ela Elb EIo 2 Lo 2 Lpa Lpb 

+ 12 Ela Elb EIo 2 Lo Lpa 2 Lpb + 12 Ela Elb EIo 2 Lo Lpa Lpb 2 

+ 4 Ela Elb EIo 2 Lpa 3 Lpb + 6 Ela Elb EIo 2 Lpa 2 Lpb 2 + 4 Ela Elb EIo 2 Lpa Lpb 3 

+ Elb 2 EIo 2 Lpa 4 

EC3 = Elb EIo Lpa 3 
EC4 = Ela EIo Lpb 3 
EC5 = Ela EIB Lo 3 
EC6 = 3 Elb EIo Lpa 2 Lpb 
EC7 = 3 Ela EIo Lpa Lpb 2 
EC8 = 3 Elb EIo Lo Lpa 2 
EC9 = 3 Ela EIo Lo Lpb 2 
EC10 = 3 Ela Elb Lo 2 Lpb 
ECU = 3 Ela Elb Lo 2 Lpa 

Tomando el mismo edificio de la figura 2 se obtuvo la Curva de Capacidad utilizando para 
cada seccion una curva de modelo completo a flexocompresion para columnas y solo a flexion 
para vigas con el metodo de Giberson, y este resultado se lo comparo con curvas obtenidas con el 
codigo ASCE 41, con el metodo de Giberson y el de la matriz K1 condensada. A continuacion se 
presenta el edificio analizado en 2 dimensiones, comparando los resultados: 



Curva V-Deformacion tope 



25C 



200 



- Con Momento Curvatura sin carqa axial y Modelo de Giberson 

- Con Momento Curvatura con Carga axial y Modelo de Giberson 
Con Momento Rotation ASCE41 y Modelo de Giberson 

■ Con Momento Rotation ASCE41 y Modelo K1 condensado 




Figura 12 Resultado de Curva de dse Capacidad 
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7.4 Matriz de rigidez condensada en coordenadas reducidas incluyendo 
deformacion por corte 

De igual manera como se vio en la seccion 7.3 podemos trabajar incluyendo la 
deformacion por corte, y con un proceso parecido a ese se puede obtener una matriz K1 , pero 
debido a su gran tamano de aproximadamente 200 terminos se expondra el proceso como 
programa. 

Los coeficientes por deformacion por corte son como se los muestra en la matriz 
acontinuacion: 



kv = (2 + beta) /{I + 2* beta); 
av = (1 — beta) /(l + 2 * beta); 

beta = 6 * EI/(A/f * G * L A 2), G = Ey/(2 * (1 + mm)), miu = .2, f = 
factor de forma de la seccion, A = Area de la seccion, y L = es la longitud del miembro: 



Programacion en Matlab: 

function[K1_condensada]=K1_eqcorte2(Ela,Elo,Elb,EAa,EAo,EAb,Lpa,Lo,Lpb,beta_a,beta_o, 
beta_b) 

%Ela,Elo,Elb: Rigidez a flexion en extremo a, en centro y en extremo b 
%EA,EAo,EAb: Rigidez axial en extremo a, en centro y en extremo b 
%Lpa,Lo,Lpb: Longitud en extremo a, en centro y en extremo b 
%beta_a,beta_o,beta_b:factor de deformacion por corte beta=6*Ely/(A/f*G*L A 2) 
%en extremo a, en centro y en extremo 

kva=(2+beta_a)/(1 +2*beta_a); 
kvo=(2+beta_o)/(1 +2*beta_o); 
kvb=(2+beta_b)/(1 +2*beta_b); 
ava=(1 -beta_a)/(1 +2*beta_a); 
avo=(1 -beta_o)/(1 +2*beta_o); 
avb=(1 -beta_b)/(1 +2*beta_b); 

sigma=(2*Elo*avo)/Lo A 2+(2*Elo*kvo)/Lo A 2; 



'2 * El/lon * kv 2 * El/lon * av 
Kl = 2 * El/lon * av 2 * El/lon * kv 
0 0 



0 
0 

EA/lon 
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K e = 

EAo EAa EAo 
— + —,0,0, -—,0,0; 
Lo Lpa Lo 

2^/2* Ela * ava 1* Ela * kva\ 

2 * sigma \ Lpa 2 Lpa 2 / 2 * EIo * avo 2 * Ela * ava 2 * EIo * kvo 

Lo Lpa ' Lo 2 Lpa 2 Lo 2 

2 * Ela * kva 2 * sigma 
Lp~a^' l^-^ma 

2 * EIo * avo 2 * Ela * ava 2 * EIo * kvo 2 * Ela * kva 2 * EIo * kvo 

Lo 2 Lpa 2 Lo 2 Lpa 2 ' Lo 

2 * Ela * kva 2 * EIo * avo 2 * EIo * kvo 2 * EIo * avo 
+ ; ,0, 



Lpa Lo 2 Lo 2 Lo 

EAo EAo EAb 
■—,0,0, — + —,0,0; 
Lo Lo Lpb 

(2* Elb * avb 2* Elb * kvb\ 
2* sigma 2* EIo* avo 2 * EIo * kvo 2 * sigma *\ Lpb 2 Lpb 2 / 2* Elb* avb 

°' To ' To 2 Lo 2 ' 0 ' To + Lpb ' Lpb 2 

2 * EIo * avo 2 * EIo * kvo 2 * Elb * kvb 

■ + - 



Lo 2 Lo 2 Lpb 2 

2 * EIo * avo 2 * Elb * avb 2 * EIo * avo 2 * EIo * kvo 2 * Elb * kvb 2 * EIo * kvo 
0 , sigma , _ , 0 ,_ _ _ — + _ _ , 

2 * Elb * kvb 
+ - 



Lpb 



Bnt=[ 0, 0, 0 

-(2*Ela*ava)/Lpa A 2-(2*Ela*kva)/Lpa A 2, 0, 0 

(2*Ela*ava)/Lpa, 0, 0 

0, 0, -EAb/Lpb 

0, (2*Elb*avb)/Lpb A 2+(2*Elb*kvb)/Lpb A 2, 0 
0, (2*Elb*avb)/Lpb, 0]; 



Kn =[ (2*Ela*kva)/Lpa, 0, 0; 

0, (2*Elb*kvb)/Lpb, 0; 

0, 0, EAb/Lpb]; 

K1_condensada=Kn-Bnt'*Ke A -1 *Bnt; 
%fin 

Tomando el mismo edificio de la figura 2, se obtuvo la Curva de Capacidad utilizando 
para cada seccion una curva del modelo del codigo ASCE 41 obtenidas con el metodo de 
Giberson y matriz k1 condensada con y sin deformaciones a corte. A continuacion se presenta el 
edificio analizado en 2 dimensiones comparando los resultados: 
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Curva V-Deformacion tope 

















































































































Con Momento Rotacion ASCE41 y Modelo de Giberson sin def a corte 

Con Momento Rotacion ASCE41 y Modelo K1 condensado sin def a corte 

Con Momento Rotacion ASCE41 y Modelo K1 condensado con def a corte 
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Figura 13 Resultados de Curva de Capacidad 
7.3 Modelo de plasticidad extendida con rigidez constante escalonada 




Figura 14 Modelo de viga de plasticidad extendida 

A traves del diagrama de masas se puede obtener la matriz de flexibilidad para el caso de 
viga presentado en la figura 14 con rigideces constantes escalonadas y la demostracion para 
obtener la matriz de flexibilidad esta detallada en la referenda Aguiar 2003. 

Para el algoritmo de calculo se detallan los siguientes terminos en un diagrama de masas 

elasticas: 



,X2 



X3 
X4 



X5 



X6 



X7 



X8 



Y' 



A10 



Figura 15 Diagrama de 
masas elasticas. 



n — 4; Numero de Escalones 
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A3 

A4 = ■ 

A5 ■ 
A6 ■■ 



1 



1) 



G/a EIo) 



Ela 



- (i - 1) 



Ela ^ 



1) 



Ela 



- (i - 1) 



n 








(eTo. ~ 


ETo) 


n 








(eTcl~ 


ETo) 
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{eTcl~ 


EIo) 


n 



XI 

;X2 
;X3 
;X4 



.Lpa . 

i ; i = 1 

n 

Lpa 

i ; i = 2 

n 

.Lpa . ? 

i ; i = 3 

n 

Lpa 
i — = 4 
n 



EIo 



;X5 ■ 



■■ — -0-1) 

Elb v J 



Lpa + Lo 

(i7 



" E,o) -,X6 = Lpa + Lo-i^;i= 1 
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A7 ■■ 



Elb 



A8 = 



A9 



1 

Elb 

1 
Elb 



(i-iy 

O-i) 
(i - 1) 



(l75 EIo) 



n 

(etb~ 


EIo) 


n 








{~EIb~ 


Elo) 



Lpb 

; XI = Lpa + Lo — i ; i = 2 

n 



; X8 = Lpa + Lo — i ; i = 3 

n 



Lpb 
n 
Lpb 



;X9 - Lpa + Lo - i ; i = 4 

n 



W ■ 



Al * XI + V 4j * (X; - #,-_!) 



;'=2 

, Al*(L-Xlf 2 v 1 ^f*((i-^-i) 3 -(i-^;) 3 ) 
^ = 3^ + Z — 3^ " 

J'=2 

De estos terminos se puede obtener la matriz de flexibilidad: 



F = 



Vy/L 



2 1/ L 



7y/L 2 



Tomando el mismo edificio de la figura 3 se obtuvo la Curva de Capacidad utilizando 
para cada seccion una curva obtenida por el codigo ASCE 41 con el metodo de Giberson, matriz 
k1 condensada con y sin deformaciones a corte y el modelo de plasticidad extendida escalonada. 
A continuacion se presenta el edificio analizado en 2 dimensiones comparando los resultados: 



Cur\a V-Deformacion tope 
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Con Momento Rotacion ASCE41 y Modelo de Giberson sin def a corte 

Con Momento Rotacion ASCE41 y Modelo K1 condensado con def a corte 

1 Con Momento Rotacion ASCE41 y Modelo de plast escalonada sin def a corte 
Con Momento Rotacion ASCE41 y Modelo K1 condensado sin def a corte 
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Figura 16 Curva Generalizada Fuerza - Deformacion de componentes o elementos 
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8. RESISTENCIA RESIDUAL DEL DIAGRAMA MOMENTO ROTACION 



Q_ 




BorA 

Figura 17 Resultados de curva de capacidad 

La resistencia residual de una seccion es la correspondiente al tramo entre el punto D y E 
marcada en el cuadro de la figura 18, esta se da debido a que, una seccion luego de perder su 
resistencia a cargas laterales en la estructura, aun tiene una reserva de resistencia para cargas 
gravitacionales, la cual valdra c*My. 

Para el calculo de una seccion que este en este tramo de la curva, se parte del sistema 
de coordenadas reducidas K1 de cada elemento, ver figura 19, se le multiplica por matrices de 
transformacion T12 para llevar al elemento a coordenadas locales K2=T12T*K1*T12 y a esta por 
las matrices T23 para llevar al elemento a coordenadas generalizadas K3=T23T*K1 *T23, 
finalmente luego de ensamblar la matriz de rigidez K de la estructura con las matrices K3 de los 
elementos, obtenemos todos los desplazamientos de los grados de libertad q=Q\K; luego las 
fuerzas en cada grado de libertad son la suma de las fuerzas de empotramiento perfecto (vector 
Q2) por cargas en el elemento y la fuerza PL = T23 * PG donde, PG = K3 * pq, y pq son los 
desplazamientos de cada grado de libertad correspondiente a este elemento, entonces las fuerzas 
seran F = Q2 + PL. 

Invirtiendo este proceso tenemos: 
F= Q2 + T23*PG 



Si: 723 = 



cosen 


sen 


0 


0 


0 


0 


—sen 


cosen 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


cosen 


sen 


0 


0 


0 


0 


—sen 


cosen 


0 


- 0 


0 


0 


0 


0 


1- 



T23 corresponde a una matriz de transformacion con grados de libertad del nudo de 
extremo numerados: primero el grado horizontal, segundo el grado vertical y tercero el giro, y los 
siguientes tres grados de libertad corresponden al extremo en el otro nudo. "sen" y "cosen" son el 
seno y coseno del angulo de orientacion del elemento con respecto a la coordenada global 
horizontal. 

Si queremos tomar en cuenta solo los momentos, tomaremos solo la columna 3 o 6 de la 
matriz T23, propio a giros en el extremo inicial o final correspondientemente, y sabiendo que F es 
c*My tenemos: 

c * My = Q2(fila 3 o 6) + K3(columna 3 o 6) * pq(fila 3 o 6) (14) 
Si la matriz de elemento en coordenadas reducidas es: 
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Conociendo ya la matriz K3 que va a estar en funcion de K1 , (K2=T12T*K1*T12, 
K3=T23T*K1 *T23) podemos despejar de la Ec.14: 

Nudo inicial pasol) Para, cuando el nudo inicial de la seccion esta entre el punto D y E del 
diagrama Momento - Rotacion: 

K 1_11 = —(Q2(fila 3) * Ion — c * My * Ion + K\_Y2 * cosen * pq(fila 2) — Kl_12 * cosen * 
pq(fila 5) + Kl_12 * Ion * pq(fila 6) — Kl_12 * pq(fila 1) * sen + K\_Y2 * pq(fila 4) * sen)/(cosen * 
pq(fila 2) — cosen * pq(fila 5) + Ion * pq(fila 3) — pq(fila 1) * sen + pq(fila 4) * sen) (15) 



Nudo final pasol) Para, cuando el nudo final de la seccion esta entre el punto D y E del 
diagrama Momento - Rotacion: 

A'l_22 = —(Q2(fila 6) * Ion — c * My * Ion + K\_Y2 * cosen * pq(fila 2) — Kl_12 * cosen * 
pq(fila 5) + Kl_12 * Ion * pq(fila 3) — Kl_12 * pq(fila 1) * sen + K\_Y2 * pq(fila 4) * sen) I [cosen * 
pq(fila 2) — cosen * pq(fila 5) + Ion * pq(fila 6) — pq(fila 1) * sen + pq(fila 4) * sen) (16) 

Nudo inicial paso2) Si se utiliza el modelo de Giberson (ver seccion 7) se puede despejar Ela 
en funcion de Elb y Elo desde la Ec.1 5 (Ela, Elb son rigideces en extremos de miembros, ver figura 
10): 

Ela = -(Elb * Lpb A 2 * Q23 + Elo * Lpb A 2 * Q23 - Elb * Lpb A 2 * cMy - Elo * Lpb A 2 * cMy)/(12 * Elb 
*p2 - 12* Elb *p5 + 6*EIo*p2 - 6*EIo*p5 + Lpb A 2 * Q23 - Lpb A 2*cMy 
+ 6 * Elb * Lpb *p3 + 6 * Elb * Lpb * p6 + 6 * Elo * Lpb * p3) 
Nudo final paso 2) O se puede despejar Elb en funcion de Ela y Elo desde la Ec.1 6: 

Elb = -(Ela * Lpb A 2 * Q26 + Elo * Lpb A 2 * Q26 - Ela * Lpb A 2 * cMy - Elo * Lpb A 2 * cMy)/(12 * Ela 
*p2 — 12* Ela * p5 + 6*EIo*p2 — 6*EIo*p5 + Lpb A 2 * Q26 — Lpb A 2*cMy 
+ 6 * Ela * Lpb *p3 + 6 * Ela * Lpb * p6 + 6 * Elo * Lpb * p6) 
Para otros modelos mas complejos que el de Giberson despejando Ela o Elb resulta en 

ecuaciones demasiado complejas de cientos de terminos, entonces podemos trabajar la matriz K1 

de la siguiente manera: 
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Nudo inicial paso 3) Para, cuando el nudo inicial de la seccion esta entre el punto D y E del 
diagrama Momento - Rotacion: 

Ya que se vio en secciones anteriores que el metodo del modelo de Giberson es un 
metodo muy aceptable, se parte de este modelo para obtener los terminos de la matriz K1 : 



Kl = 



Kill Kill 0 
Kill Kl_22 0 
0 0 Kl_33 



Del Modelo de Giberson (Seccion 7.2): 

Kl_12 = (6 * EI a * Elb) /(Ion * (EI a + Elb + EIo)) 

Kl_22 = (6 * Elb * (Ela + EIo))/(lon * (Ela + Elb + EIo)) 

Si Ela=q*EI, Elb=p*EI y Elo=EI 

6 * p * q * EI 



Kl 12 = 



Kl 22 = 



I * (p + q + 1) 

6 * p * EI * (q + 1) 6 * p * q * EI 



(q + 1) 



Kl 12 = 



L * (p + q + 1) 
Por lo tanto: 



L * (p + q + 1) 



(c/ + D 



Kl 22 



Tomando esta matriz K1 como valida para todo modelo, no solo el de Giberson, e 
invirtiendo el proceso con Kl_22 = Elb * (2 + beta)/(l + 2 * beta), se tiene: 



Kl = 



(Ver Ec. 15) 



-Kl 22 



Kl 22 



0 



(q + D 

6* Elb* EIo (1 2 + beta \ 
L * (Elb + EIo) * \2 * 1 + 2* beta) 

0 Rigidez Axial 



En la resistencia residual que se tendra en la seccion, para alcanzar un momento igual a 
c (0.2*My) se tendna que usar una rigidez tendiendo a ser cero, q«0. 



Kl = 



(Ver Ec. 15) 
0 
0 



6*EIb*EIo 
L*(EIb+EIo) 



0 

■G- 

0 



2+beta \ 
l+2*betaj 



0 
0 

Rigidez Axial 



(Ec.17) 



beta es el coeficiente que considera deformacion por corte ver seccion 7.4, beta = 6 * 
EI/(A/f*G*L A 2). 

Nudo final paso 3) Para, cuando el nudo final de la seccion esta entre el punto D y E del 
diagrama Momento - Rotacion: 
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' b*Ela*Elo (1 2+beta \ 


0 


0 












L*{EIa+EIo) \2 l + 2*beta) 


(Ec.18) 


Kl = 


0 


(Ver Ec. 16) 


0 




0 


0 


Rig idez Axial 





9. RELACION APROXIMADA ENTRE MOMENTO CURVATURA Y 
MOMENTO ROTACION 



9.1 Metodo longitud de la rotula plasties 




En este metodo la matriz de rigidez de una columna es 
formulada por la combinacion de un elemento de columna elastica y un 
elemento de resorte rotacional conectado en cada extremo, como se 
ve en la figura 19. 

Los resortes rotacionales y la rigidez elastica a flexion de la 
columna se comportan con dos resortes en paralelo. La rigidez del 
resorte rotacional es gobernado por la curva momento - rotacion de 
una rotula con longitud Lp. 

La matriz de rigidez del miembro puede ser derivada usando la 
teoria modificada de la pendiente - deflexion: 

AM £ = - [4(A0i - A9i H ) + 2(A0y - A0;' H )] (19) 



Donde: 



Figura 19 

Combinacion de columna 
elastica y resortes 
rotacionales 



AM) = — [4(A0i - A8i H ~) + 4(A0) - A9j H )] 



AO- AMi A/T AM J 

A9i H = — y A0t w = — - 

" Ri ' 3 " Rj 



(20) 



(21) 



A9i H y A8i H son las rotaciones plasticas en los extremos del elemento de la figura 20. Ri y 
Rj son la rigidez inelastica de la articulacion plastica en los extremos del elemento de la figura 20, y 
pueden ser obtenidos desde las pendientes de la curva M - 8 H correspondientes a los extremos del 
elemento de la figura 20. Resolviendo Ec.1 9 y Ec.20 para Mi y Mj se obtiene: 

AMi = - [Sii *A9i + Sij* A9j] (22) 



AM; = j [Sij * A6i + Sjj * A0j] 



(23) 
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Donde, 



fAMij _ E/ [Sit Sin (A6i) 
\AMj\-T[sij SjjlUej] 



Sii = [4 + ^]/R* 
Sij = £ 

V LRU V Lfi;/ Vl/ Vflifl;/ 



(28) 



La rotacion total en el tope de la rotula plastica de la columna puede ser expresada como 
0 = 0 y + 0 H , donde 0 y es la rotacion de fluencia en el tope de la rotula y 0 H es la rotacion plastica 
en el tope de la rotula. Basado en la relacion M - <j), 0 y y 0 H pueden ser calculadas como sigue: 

rip r ip M f Lp My/L-y\ My\ r Lp 1 r Lp 

* dy =h K dy= l k(— ) dy= w\l dy -zJ„ ydy 

e„ = <$> H Lp 
o H = (4>- <$> y )Lp 



Lp 



(30) 



Donde cj) es la curvatura total correspondiente al momento, en el fondo de la columna, 
asumiendo que la curvatura plastica, (cf> - c)) n ), es uniformemente distribuida a lo largo de la 
longitud plastica, Lp. 



_ (fl-fly) My 
^ Lp Ele 



(31) 



Para el diagrama de Momento - Rotacion generalizado del codigo ASCE 41 , figura 21 , la 
curvatura en el punto C sera: 



4>c = <K ; 



({By + a) - 6 y ) My 



Lp 



Ele 



(32) 
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Figura 20 Curva generalizada Fuerza-Deformacion de componentes o elementos 

Por geometna del diagrama Momento-Curvatura figura 21 podemos determinar: Mu = 
(Ou - <t>y) * EIp + My, y reemplazando con la Ec. 35, tendremos: 




0u 0 



Figura 21 Diagrama Momento-Curvatura Bilineal 

Mu = (c() c - c()y) * EIp + My 



Mu 



= fe + *y)-*yJ*Eip + My 



Mu = — * EIp + My 



(33) 



V-Desplazamiento en tope 











































































































Modelo teta H =fi H *Lp 

Modelo Directo Momento-Rotacion de ASCE 41 
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Figura 22 Curva de capacidad con modelo 8 H = § H Lp y con Modelo de Momento Rotacion directo. 
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10. 



CONCLUSIONES 



• Existen diversos modelos, simples y complejos, lo que resulta en tiempo de calculo 
computacional. Teoricamente el modelo visto en la seccion 3 da resultados mas reales 
pero conlleva un gran esfuerzo computacional, se recomienda utilizar un modelo mas 
simple como el modelo Momento rotacion del codigo ASCE 41, aunque, se dan resultados 
mas conservadores (ver figura 13), el codigo da parametros sacados no solo de teoria de 
materiales sino tambien de experimentacion, ademas se tiene el respaldo de un codigo tan 
importante como el ASCE 41, se recomienda usar estas curvas con una pendiente de 
fluencia del 3% de la pendiente elastica como especifica este codigo. 

• Enfocandonos en el tipo de analisis para la matriz de rigidez en coordenadas locales de un 
elemento todos los estudiados aquf son muy aproximados en sus resultados. Aunque sea 
cierto un cambio no lineal de rigidez a flexion El en las diferentes secciones de un 
elemento, desde una seccion en rango inelastico hasta una seccion en rango elastico (ver 
figura 8), usar un modelo real de cambios de rigideces en la rotula plastica como la 
mostrada en la figura 25, sera equivalente a usar un modelo como el de la figura 26 que es 
el usado en el modelo de la matriz de elemento K1 condensada. Este modelo mencionado 
a demas tiene la ventaja de poder usar rigideces axiales variables a lo largo de un 
elemento. 



M 




M 



Figura 23 Variacion de rigidez real en 
longitud de rotula plastica de un elemento 



Figura 24 Variacion de rigidez de modelo 
bilineal en rotula plastica 



Es importante conocer todos estos procesos para la obtencion de curvas de capacidad, y 
usar con seguridad un modelo para los analisis, sabiendo los defectos, virtudes y ventajas que se 
tendra en el resultado final. 
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